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SUMMARY 


A  conventional  solution  of  the  Helmholtz  or  time-reduced  wave  equation  is 
a  simple  product  of  functions  that  contain  one  coordinate  variable  in  each. 
An  unbounded  set  of  solutions  that  are  not  separable  into  simple  products 
of  single -variable  functions  has  been  partially  examined  for  applicability 
to  vibrational  problems.  Applications  to  scalar  usage  have  been  found, 
and  illustrations  including  shapes  and  frequencies  for  membranes  and  an 
acoustic  cavity  are  reported.  Efforts  to  make  application  to  vector  usage 
are  described,  as  are  numerous  mathematical  properties  that  have  been  dis¬ 
covered  in  the  course  of  the  'work.  It  is  concluded  that  vibration  on  or 
within  some  new  shapes  can  now  be  calculated  exactly  with  functions  formed 
of  the  nonseparable  solutions  added  to  separable  solutions.  It  is  also 
concluded  that  simplifications  in  the  mathematics  and  additional  appli¬ 
cations  await  the  effort* 
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FOREWORD 


This  is  a  comprehensive  report  of  an  investigation  which  has  led  more 
deeply  into  mathematics  than  was  originally  planned.  In  order  that  the 
reader  who  has  not  the  time  to  dwell  upon  the  mathematical  side  may  ob¬ 
tain  a  qualitative  description  of  the  work,  the  report  has  been  organized 
in  such  a  way  as  to  give  comprehensive  coverage  without  formulas  in  the 
Summary,  the  Introduction,  the  Course  of  the  Work  Covered  by  Contract, 
the  Conclusions  and  the  Recommendations.  Fundamental  mathematics  that 
is  the  essential  groundwork  for  the  positive  findings  is  then  given  and 
it  leads  into  the  sections  which  give  details  of  the  applications.  Other 
properties  of  the  nonseparable  solutions  were  discovered  which  do  not  lead 
directly  to  the  reported  applications,  and  these  are  detailed  in  Appendixes 
I  through  VII. 

Principal  personnel  in  this  investigation  have  been  Dr.  D.  S.  Moseley  as 
Project  Leader  and  James  M.  Render.  Mr.  Render  calculated  the  patterns, 
led  in  the  vector  effort  and  has  contributed  importantly  to  the  present 
level  of  understanding  and  to  the  preparation  and  editing  of  this  report. 
This  report  was  typed  by  Mrs.  G.  J.  Fowler. 

The  work  period  covered  by  contract  was  28  June  1965  to  28  April  1966. 

The  contract  was  DA  UU-177-AMC-3i*2(T)  between  the  U.  S.  Army  Aviation 
Materiel  Laboratories,  Fort  Eustis,  Virginia,  and  Vitro  Laboratories, 

Silver  Spring,  Maryland.  Title  of  the  contract  was  "Application  of 
Nonseparable  Wave  Equation  Solutions". 
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INTRODUCTION 


VIBRATION  AS  WAVES 


Vibration  is  oscillatory  motion  of  matter  in  response  to  a  transient  or  a 
periodic  disturbance.  Oftentimes  the  disturbance  that  is  perceived  by  a 
particular  portion  of  matter  has  originated  at  considerable  distance  and 
has  been  propagated  to  that  point  as  a  wave.  In  these  cases,  provided  am¬ 
plitudes  are  sufficiently  anall,  the  propagation  through  the  connecting 
medium  is  described  by  the  wave  equation  written  for  fluids  and  by  the 
wave  and  flexural  wave  equations  written  for  solids.  These  are  partial 
differential  equations,  one  or  both  of  which  permit  prediction  of  depend¬ 
ent  variables  such  as  pressure,  stress,  stress  moment,  strain,  displace¬ 
ment,  particle  velocity,  and  acceleration  anywhere  in  the  medium  and  at 
any  moment  of  time. 

For  ease  in  handling  these  equations,  one  usually  introduces  the  assumption 
that  the  dependent  variable  is  a  sinusoidal  function  of  time.  This  changes 
wave  and  flexural  wave  equations  into  second  and  fourth  order  partial  dif¬ 
ferential  equations,  respectively,  in  position  coordinates  only,  with  fre¬ 
quency  appearing  as  an  arbitrary  parameter.  The  wave  equation  so  reduced 
is  called  the  Helmholtz  wave  equation. 

SEPARABLE  AND  NONSEPARABLE  SOLUTIONS 


Functions  which  are  solutions  of  the  Helmholtz  wave  equation  have  been 
known  for  a  century  and  more.  They  are  derived  by  separation  of  vari¬ 
ables,  a  method  taught  to  engineering  and  science  students  in  colleges 
everywhere.  The  method  assumes  that  the  solution  is  a  function  that  is 
itself  a  simple  product  of  functions  each  containing  one  independent  vari¬ 
able.  By  this  method  the  derivation  is  converted  from  one  of  solving  a 
partial  differential  equation  in  two  or  three  dimensions  to  one  of  solving 
two  or  three  ordinary  differential  equations.  In  view  of  the  nature  of 
this  solution  and  to  make  a  distinction  between  it  and  those  to  follow, 
we  have  called  it  the  separable  solution  of  the  Helmholtz  wave  equation. 

The  nonseparable  solutions  of  the  Helmholtz  wave  equation  with  which  this 
report  is  concerned  are  mathematical  functions  which  have  appeared  in  the 
past  U  years.  They  arose  during  the  course  of  an  attempt  to  formulate 
in  closed  form  the  natural  modes  of  an  isotropic  right  circular  solid 
cylinder.  No  useful  role  for  them  was  found  during  that  work,  but  they 
were  later  gathered  up  and  submitted  for  publication  (reference  ii). 

Preparation  was  then  made  for  the  effort  that  is  being  reported  here. 

PREPARATION  FOR  THIS  CONTRACT 


It  was  reasoned  that  new  solutions  of  an  equation  that  is  used  in  mechanics, 
acoustics,  radio,  microwaves,  and  light  would  have  utility.  It  was  postu¬ 
lated  that  these  solutions  in  mechanics,  for  example,  might  legitimatize 
experimentally  observed  modes  which  theoreticians  dismiss  as  originating 
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in  poorly  controlled  homogeneity  of  medium  and  poorly  known  boundary  con¬ 
ditions. 

It  was  decided  to  change  from  circular  cylindrical  coordinates  of  the  pub¬ 
lished  paper  to  rectangular  coordinates.  The  aptness  of  the  two-dimensional 
nonseparables  to  membrane  vibration  was  then  apparent,  since  the  motion 
being  predicted  was  in  a  single  direction,  namely,  normal  to  the  plane  of 
the  membrane,  and  the  dependent  variables  were  the  minimum  number,  namely, 
two,  that  were  required  for  the  existence  of  the  nonseparable  solution. 

Nodal  patterns  for  two  of  the  new  functions  were  quickly  generated,  and 
certain  features  of  these  patterns  were  noted  as  being  identical  with 
certain  features  of  Chladni  figures  republished  by  VJaller  (reference  7). 
Chladni  figures  are  sand  patterns  formed  on  vibrating,  free-edged  plates. 


Figure  1.  Selected  Chladni  Figures  on  Square  Plates  (reference  7). 

A  program  was  then  outlined  that  would  successively  investigate  appli¬ 
cations  in  four  models,  namely,  membranes,  acoustic  cavities,  electro¬ 
magnetic  cavities,  and  elastic  solids.  These  were  so  ordered  because  it 
was  felt  that  this  was  their  order  of  difficulty  and  because  it  was 
thought  unwise  to  undertake  only  the  first  and  last  without  the  experi¬ 
ence  to  be  gained  in  executing  the  intermediate  ones.  First  attention 
in  each  was  to  be  given  to  finding  if  new  modes  for  shapes  already  known 
would  be  given  by  the  nonseparables.  For  each  model  there  were  guides 
proposed  as  to  number  of  frequencies,  number  of  regular  shapes,  highest 
order  of  solution  to  be  considered,  etc. 

Between  proposal  and  contract,  a  factor  of  considerable  importance  was 
discovered.  There  was  found  a  generative  operator*  by  which  the  simplest 
nonseparable  could  be  obtained  from  the  known  separable,  and  this  operator 
proved  capable  of  generating  ever  higher  orders  of  nonseparable  through 
successive  application. 


^Credit  for  recognizing  this  operator  belongs  to  Dr.  B.  R.  Levy,  Mathematics 
Branch,  Office  of  Naval  Research,  Washington,  D.  C. 
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COURSE  OF  THE  WORK  COVERED  BY  CONTRACT 


THE  SEARCH  FOR  APPLICATIONS 


Under  Contract  DA  bli-177-AMC-3U2(T)  entitled  ’’Applications  of  Nonseparable 
Wave  Equation  Solutions”,  the  U.  S.  Army  Aviation  Materiel  Laboratories  has 
been  supporting  a  research  effort  that  began  as  sr.  effort  to  discover  if 
applications  exist.  The  initial  model  was  a  membrane. 

It  was  soon  apparent  that  nonseparable  solutions  for  membranes  were  not 
producing  new  modes  in  old  shapes,  for  calculated  nodal  lines  on  membranes 
bore  no  resemblance  to  regularly  bounded  patterns.  Therefore, that  goal 
was  shelved  in  favor  of  achieving  a  nodal  pattern  of  any  shape  that  would 
demonstrate  one  or  more  closed  areas.  Nodal  lines  that  do  not  close  on  an 
infinite  membrane  are  interesting  but  not  practical. 

When  the  schedule  called  for  it,  the  program  moved  from  membranes  to 
acoustical  cavities,  that  is,  from  two  to  three  dimensions.  More  patterns 
were  evaluated.  At  length  an  x-y  pattern  in  a  z-plane  did  demonstrate 
closed  areas,  and  study  of  the  solution  which  had  produced  them  revealed 
the  secret:  in  effect  the  plot  was  of  a  function  very  similar  to  a  mem¬ 
brane  function  compounded  of  a  nonseparable  added  to  a  strong  proportion 
of  separable.  A  return  to  membranes  was  made,  and  closed  areas  were  im¬ 
mediately  demonstrated.  Applied  to  acoustical  cavities,  the  principle  of 
using  a  nonseparable  to  perturb  a  separable  solution  was  successful  in 
generating  closed  volumes. 

In  preparation  for  the  third  model,  examination  was  made  of  nonseparable 
solutions  as  scalar  potentials  from  which  vectors  could  be  derived.  These 
were  then  studied  in  relation  to  boundary  conditions  on  the  cavity  walls, 
and  this  included  a  return  to  acoustical  cavities  with  normal  particle 
velocity  to  be  zero  on  the  wall.  The  mathematical  meaning  of  two  boundary 
conditions,  electric  and  magnetic,  needing  to  be  met  on  the  same  wall  was 
investigated,  and  electric  and  magnetic  potentials  were  constructed. 

Attention  then  turned  to  elastic  solids  where  two  conditions  are  required 
on  each 'surface  or  edge.  Zxtensional  and  flexural  vibration  in  thin  plates 
was  examined. 

SUMMARY  OF  RESULTS 


Application  to  vibrating  membranes  has  been  found.  With  the  aid  of  the  new 
solutions,  frequency  and  amplitude  distribution  on  shapes  never  before  ana¬ 
lyzed  have  been  found  exactly  and  without  recourse  to  representation  by 
Fourier  series.  It  appears  that  the  exact  solution  of  any  shape  that  is  a 
continuous  distortion  of  a  square  or  rectangle  is  within  grasp. 

Application  to  an  acoustic  cavity  with  pressure-node  walls  has  been  demon¬ 
strated.  Frequency  and  amplitude  distribution  of  pressure  within  the  volume 
can  new  be  calculated.  Wall  shape  in  this  demonstration  is  a  smoothly  dis¬ 
torted  cube. 
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Application  to  an  electromagnetic  cavity  of  perfectly  conducting  walls  was 
discovered  during  the  preparation  of  this  report  and  has  not  been  well  con¬ 
firmed.  It  appears  that  nonseparables  of  first  order  can  serve  as  electric 
and  magnetic  potentials  from  which  are  obtained  electric  and  magnetic  fields 
that  predict  fields  within  a  hollow  rectangular  parallelepiped.  Dimensions 
of  the  cavity  must  be  proportional  to  any  set  of  three  integers. 

There  has  been  no  success  with  the  acoustic  cavity  under  the  condition  of 
a  rigid  wall  and  only  the  foregoing  success  with  the  electromagnetic  cavity. 
These  models  have  in  common  the  boundary  conditions  that  tangential  and/or 
normal  components  of  vector  field  quantities  are  specified,  and  much  effort 
has  been  expended  in  trying  to  use  nonseparable  solutions  to  represent  suit¬ 
able  vector  fields. 

The  fourth  model  is  the  elastic  solid.  Two  modes  of  motion  in  a  thin 
rectangular  plate  were  examined  in  the  last  weeks  of  the  contract.  The 
first  was  extensional  vibration,  which  is  motion  parallel  to  the  surface 
of  the  plate.  Nonseparables  appear  unable  to  satisfy  the  pair  of  stress 
conditions  that  apply  at  each  edge.  It  may  be  noted  that  separable  s  do 
not  satisfy  them  either .  The  second  was  flexural  vibration,  or  motion 
normal  to  the  plane  of  the  plate.  Thir  examination  was  cut  short  by 
schedule,  but  it  appears  that  nonseparables  in  trigonometric  and  hyperbolic 
forms  can  be  taken  together  to  satisfy  boundary  conditions  along  a  plate 
edge.  The  prediction  is  that  the  edge  shape  will  in  general  be  a  distorted 
rectangle,  just  as  the  edge  shape  of  the  membrane  has  been  found  to  be. 
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MATHEMATICAL  FOUNDATIONS 


THE  HSIMHOLTZ  WAVE  EQUATION 


For  amplitudes  sufficiently  small,  the  vector  differential  equation  for 
wave  motion  is 


VSA 


1.^1 

ca  at3 


(i) 


where  the  vector  A  is  the  dependent  variable,  7  is  the  vector  differential 
operator,  c  is  the  phase  propagation  velocity,  and  t  is  time.  The  7  op¬ 
erator  taken  twice  is  the  Laplacian,  which  in  rectangular  coordinates  is 


The  velocity  c  is  that  velocity  which  is  appropriate  to  wave  type  and 
medium,  examples  being  velocities  of  sound  pressure  in  a  gas,  dilatational 
or  shear  stress  and  strain  waves  in  a  solid,  and  electromagnetic  wave 
components  in  a  vacuum. 

The  Helmholtz  wav.j  equation  follows  immediately  from  Equation  (l)  if  A  is 
assumed  to  be  a  vector  function  which  varies  sinusoidally  in  time.  Thus, 
let 


4  T  «iwfc 
A  ■  A  e 

o 


(2) 


in  which  Aq  is  the  amplitude  of  vector  1,  i  is  and  u>  is  the  angular 
frequency  of  the  time  variation.  Then  Equation  (1)  becomes 


(*s  ♦  k“)Ao  -  0  ,  (3) 

in  which  k  •  u>/c. 

Since  the  amplitude  Aq  is  a  vector,  in  component  form  it  is 

A'o  ■  Ux  *  JAy  ♦  Uz  • 


and  upon  substitution  in  Equation  (3)  there  arise  three  independent  scalar 
Helmholtz  equations  which  are 

(V*  ♦  kS)Ax  *  0 

(VS  ♦  k“)Ay  -  0 

(?a  ♦  k‘)kz  -  0  (li) 
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SEPARABLE  SOLUTIONS 


Separable  solutions  of  scalar  Helmholtz  wave  equation  are  those  which  con¬ 
sist  of  single,  simple  products  of  functions  containing  one  variable  each* 
For  example,  the  separable  solutions  of  the  last  of  Equation  (h)  are,  in 
three  dimensions. 


A  ■ 

z 


(5) 


where 


•  9  2  ,3  ,2  /  ✓  \ 

k  «  a  ♦  b  +  d  .  (6) 

Each  brace  contains  two  functions  which  may  be  summed  in  any  proportion  but 
which  depend  upon  one  coordinate  variable  only.  The  function  A  is  the 
product  of  the  three  braces. 

An  alternative  way  of  writing  Equation  (5)  is  with  the  use  of  arbitrary 
phase  constants,  as 

A  -  sin  (ax  +  a)  sin  (by  +  3)  sin  (dz  +  6)  .  (7) 

z 

NONSEPARABLE  SOLUTIONS  IN  TWO  DIMENSIONS 

Let  us  consider  the  scalar  Helmholtz  wave  equation 

(Va  +  k3)W-0,  (8) 

in  which  W  is  a  function  of  x  and  y  and  not  of  z. 

This  will  be  recognized  as  the  wave  equation  for  motion  on  a  membrane  when 
W  is  the  amplitude  of  transverse  displacement  (or  velocity  or  acceleration) 
everywhere  on  the  membrane. 

The  separable  solution  of  Equation  (8)  is 

-  sin  (ax  +  or)  sin  (by  +  B),  (9) 

for  which  the  frequency  equation  is 

ka  -  a3  ♦  ba.  (10) 

In  addition  to  the  separable  solution,  there  is  a  set,  infinite  in  number, 
of  nonseparable  solutions  of  Equation  (8),  for  each  of  which  the  frequency 
equation  is  the  same  as  Equation  (10).  Each  nonseparable  contains  the  co¬ 
ordinate  variables  as  explicit  multiplying  factors,  and  since  the  highest 
exponent  or  power  that  is  present  on  such  factor  is  a  unique  characteristic 
of  that  solution,  it  has  been  used  to  declare  the  order  of  the  solution. 


Nonseparable  solutions  of  the  first  three  orders  in  two  dimensions  are 

*  bx  cos  (ax  +  a)  sin  (by  +  0)  -  ay  sin  (ax  +  a)  cos  (by  +  B),  (11) 

>  -  {[(bx)2  +  (ay)3]  sin  (ax  +  or)  sin  (by  +  0)  + 

by  sin  (ax  +  rv)  cos  (by  +  8)  + 
ax  cos  (ax  +  or)  sin  (by  +  B)  + 

.  o 

2abxy  cos  (ax  +  or)  cos  (by  +  B)j,  (12) 

W(3)  *  j^3(bx)a  +  (ay)2  +  l]  ay  sin  (ax  +  a)  cos  (by  +  0)  + 

2  2 

3ab  (x  -  y  )  sin  (ax  +  cr)  sin  (by  +  0)  ♦ 

2  2 

3xy  (a  -  b  )  cos  (ax  +  cr)  cos  (by  +  0)  - 

[^3 (ay)2  +  (bx)3  +  l]  bx  cos  (ax  +  cr)  sin  (by  +  B).  (13) 

A  nonseparable  solution  of  order  n  +  1  can  be  obtained  from  that  of  order  n 
by  employing  the  generative  operator  The  relationship  is 

w(n+1>-°2DW<n>,  n  *  0  (1U) 

where 

°2d  "  ^  “  *re*  (1^ 

It  is  easily  verified  that  the  separable  solution  may  be  considered  a  non¬ 
separable  solution  of  order  zero  and  hence  a  member  of  the  set.  This  is 
also  clear  in  Appendix  I  and  Appends  II,  which  present  exponential  and 
variant  forms  of  two-dimensional  solutions. 
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APPLICATION  TO  MEMBRANES 


Motion  at  right  angles  to  the  plane  of  a  membrane  is  the  simplest  two- 
dimensional  application  of  the  Helmholtz  wave  equation.  This  is  true 
because  (1)  the  motion  can  be  treated  as  a  scalar,  and  (2)  a  single  coiw 
dition  of  zero  displacement  at  the  boundary  completes  the  problem.  Sep¬ 
arable  solutions  of  Equation  (8)  in  appropriate  coordinates  have  long 
been  known  to  satisfy  this  condition  for  simple  shapes  such  as  rectangles, 
circles,  ellipses,  etc.,  yielding  a  spectrum  of  resonant  frequencies  for 
each  shape  and  the  amplitude  distribution  for  each  frequency. 

Out  of  respect  for  this  background,  nonseparable  solutions  (W^n^  in  rec¬ 
tangular  coordinates  with  n  ■  1,  2,  3)  were  first  studied  in  relation  to 
membranes.  This  has  resulted  in  three  findings,  of  which  two  are  negative 
and  one  is  positive.  The  negative  findings  are  (1)  that  nonseparable 
solutions  have  yielded  no  new  modes  for  the  rectangular  shape  that  is 
already  solved  by  the  separable  solution,  and  (2)  that  nonseparable  solu¬ 
tions  taken  singly  do  not  describe  the  vibration  of  a  membrane  of  any 
finite  shape.  Ety  the  latter  we  are  saying  that  nodal  patterns  of  non- 
separables  of  first  and  higher  order  contain  no  closed  areas.  A  col¬ 
lection  of  nodal  patterns  for  various  nonseparables  having  a  «=  b  and 
taken  singly  is  presented  in  Appendix  III. 

DISCOVERY  OF  CLOSED  AREAS 


The  positive  finding  is  that  a  nonseparable  of  second  order  can  be  added 
to  a  separable  to  yield  frequency  and  amplitude  within  a  closed  boundary 
that  is  a  smooth  distortion  of  the  boundary  given  by  the  separable  alone. 

If  T)  be  defined  as  the  oscillatory  displacement  of  any  point  of  the  mem¬ 
brane,  then  we  may  write 

7]  =  [W<°>  +  CW^2)]  cos  cut,  (16) 

where  C  is  a  constant  and  the  bracket  is  the  amplitude  of  the  displacement. 
The  amplitude  is  thus  given  for  all  values  of  x  and  y.  Zero  amplitude 
occurs  at  all  points  (x,,  y. )  which  solve  the  two-dimensional  transcendental 
equation 

w(°)  +  cw(2)  .  (17) 

Plot  of  all  points  (x^  y^)  in  x-y  space  is  the  nodal  pattern  of  the  dis¬ 
placement  Tj.  As  would  be  expected,  each  nodal  line  separates  displacements 
that  differ  from  each  other  by  180#  in  phase. 

Closed  areas  have  been  found  in  three  evaluations  of  Equation  (17)  using 

or  *  8  ■  0,  and  a  «  b  -  1,  and  C  ■  (12. 5tt2)-1  =  .0081,  C  *  (6.5tt2)-1  «=  .01^6, 

and  C  «  (2.5TT3)”1  ■  ,0U05.  These  patterns  together  with  the  one  for  C  =  0 
are  presented  as  Figure  2. 


8 


c  =  o 


Figure  2.  Nodal  Patterns  on  Infinite  Membrane  for  Sums  of  Separable 
and  Nonseparable  Solutions.  (Note  Closed  Areas.  Function 


is 


+  cv/2) 


with  a  -  b  ■  1  and  a  ■  B  ■  0.) 


9 


Several  features  of  these  patterns  may  be  noted.  The  outermost  closed  bound¬ 
aries  are  not  similar  in  the  geometric  sense  of  the  word.  The  number  of 
closed  areas  decreases  as  the  proportion  of  nonseparable  to  separable  in¬ 
creases,  until  it  is  zero  at  a  value  of  C  little  larger  than  G  »  ,0U05.  The 
distortion  of  the  square  that  is  the  building  block  for  C  =  0  increases  with 
distance  from  the  origin  in  each  of  the  lower  patterns.  The  distortion  of 
the  square  that  corners  at  the  origin  increases  as  C  increases.  A  character¬ 
istic  of  the  three  patterns  formed  by  the  addition  of  nonseparable  solution 
to  separable  solution  is  that  most  of  the  nodal  lines  of  the  pattern  do  not 
close  on  themselves  to  enclose  finite  areas  of  membrane. 

PRACTICAL  USE  OF  THESE  PATTERNS 


Any  nodal  line  in  the  patterns  of  Figure  2  can  become  the  edge  of  a  mem¬ 
brane  by  supposing  that  a  clamp  conforming  to  the  line  is  impressed  there. 
If  a  closed  contour  is  selected,  it  can  serve  as  the  boundary  cf  a  membrane 
of  finite  extent.  Displacement  at  any  instant  within  and  on  that  boundary 
is  given  by  Equation  (16),  and  frequency  is  given  by 


o> 


a ,  a  .  a  N 
c  (a  +  b  ) 


2ca, 


(18) 


for  every  pattern  since  a  ■  b  ■  1  throughout  Figure  2. 


The  following  procedure,  illustrated  by  Figure  3,  may  be  used  to  compute 
frequency  of  membrane  vibration  within  any  closed  contour  shown  in  Figure  2. 
First,  we  remove  the  requirement  that  the  propagation  constants  are  equal 
to  unity.  Second,  the  nodal  patterns  are  considered  to  be  points  (ax.^,  ay^) 

plotted  in  X-Y  space,  where  X  *  ax  and  Y  «  ay.  Third,  the  pattern  for 
C  ■  0  is  superposed  upon  the  one  containing  the  selected  pattern  by  causing 
the  axes  to  coincide.  Frequencies  of  the  two  patterns  are  identical  in 
accordance  with  Equation  (18),  and  that  of  the  square  pattern  is  the  refer¬ 
ence  given  by 


„  a  a 
2c  TT 

.3  * 


(19) 


where  L  is  the  side  of  the  square.  Fourth,  with  aL  the  X-side 
let  aS  be  the  X-side  of  the  desired  contour.  Note  that  aL  =  tt 
Fifth,  convert  axes  from  X-Y  to  x-y  and  convert  square  of  side 
of  side  S.  Propagation  constant  a  ■  tt/L  then  becomes  a/  ■  tt/S 
frequency  of  the  desired  contour  is  given  by 


of  the  square, 
and  aS  -  Mtt. 

L  into  square 
»  TT/toL.  Thus 


U) 


0  a  3. 
2c  tt  ] 


M2 


(20) 


Figure  3  has  been  prepared  to  illustrate  this  procedure.  The  upper  and 
lower  left  diagrams  are  superpositions  of  the  square  pattern  upon  the  de' 
sired  contour  in  X-Y  space.  The  upper  and  lower  right  diagrams  show  tho 
square  converted  to  side  S. 


IS-'Sfc- 
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3.4577T 


Figure  3.  Examples  of  Construction  Used  To  Compute  Frequency  of 
Contoured  Membrane, 
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Frequency  ratio  M  has  been  computed  for  all  patterns  of  Figure  2.  These 
are  displayed  in  Figure  U,  where  side  lengths  of  contours  chosen  from 
Figure  2  have  been  made  equal  to  clarify  the  presentation. 

Also  included  in  Figure  are  six  contours  which  are  predicted  for  an 
isosceles  right  triangle,  solutions  for  which  have  long  been  known  (refer¬ 
ence  page  80,  and  reference  3,  page  755)  and  require  no  nonseparable 
solutions.  They  were  chosen  for  presentation  because  of  their  resemblance 
to  patterns  generated  in  this  program.  Values  of  M  have  also  been  calcu¬ 
lated  for  them,  and  these  are  presented  beneath  each  triangle.  The  multi¬ 
ple  values  beneath  the  triangles  correspond  to  the  different  contours  in 
each  and  are  to  be  compared  as  directed  by  the  arrows. 

It  should  also  be  noted  that  the  M-ratios  for  the  in-triangle  contours  are 
greater  in  every  case  than  those  of  the  corresponding  contours  derived  from 
nonseparables.  However,  the  difference  is  perhaps  too  small  to  be  of  en¬ 
gineering  significance,  which  suggests  that  the  triangle  model  with  a  cor¬ 
rection  factor  might  be  useful  where  applicable  in  design  work. 

It  is  expected  that  patterns  within  modified  rectangles,  i.e.,  for  a  /  b, 
can  be  evaluated  with  the  aid  of  nonseparables,  whereas  the  only  triangle 
solution  of  which  we  know  corresponds  to  a  =  b. 

Note  snould  also  be  taken  of  the  fact  that  none  of  the  higher  frequencies 
found  for  patterns  of  this  contract  may  properly  be  called  overtones  since 
the  outermost  contours  are  not  similar  in  the  geometric  sense  to  those  of 
the  lowest  frequency.  How  to  make  them  similar  remains  a  challenge. 


M  =  1 .000 


C  =  0.008l  C  =  0 

M  =  3.457  M  =  l. 032,  2.081,3.536 

n  =  2 


Figure  h •  Closed  Contours  and  Their  Frequencies  in  Ratio  to  That 
of  Square  of  Same  Side  Length. 
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APPLICATION  TO  ACOUSTICAL  CAVITIES 


DISCOVERY  OF  CLOSED  VD LUKES 


The  ability  of  nonseparable  solutions  to  generate  nodal  surfaces  which  en¬ 
closed  finite  volumes  was  in  doubt  for  some  time.  It  was  at  length  demon¬ 
strated  that  the  addition  of  a  small  portion  of  nonseparable  to  a  separable 
was  at  least  one  way  to  produce  closed  nodal  surfaces,  just  as  had  been 
found  for  the  two-dimensional  case  of  the  membrane. 

Contributing  to  the  perplexity  was  the  discovery  that  in  three  dimensions 
more  than  one  nonseparable  solution  of  each  order  exists.  At  least  three 
appear  to  have  been  isolated.  Via  a  three-dimensional  generative  operator 
(see  Apoendix  IV)  discovered  early  in  the  contract,  nonseparable  solutions 
of  first  and  second  order  can  be  written.  A  symmetric  nonseparable  of 
second  order  was  synthesized  after  study  of  results  with  the  operator,  yet 
it  appears  distinct  from  its  progenitor.  The  third  was  a  nonseparable  of 
first  order  which  was  synthesized  directly  from  study  of  properties  that  a 
function  must  have  to  be  a  solution  of  the  Helmholtz  wave  equation,  and  it 
seems  unrelated  to  the  first  two. 

The  first  and  third  of  these  are  presented  in  Appendix  IV,  and  nodal 
patterns  of  the  former  are  presented  in  Appendix  V.  The  second  was  chosen 
for  addition  to  zero  order  and  is  given  here. 

The  chosen  function  is  a  symmetric  function  when  a  *  b  »  d  «  1  and  when 
or  ■  R  ■  6  b  o,  since  it  then  transforms  into  itself  when  its  variables  are 
interchanged  in  pairs.  It  is 

■  (xa  +  ya  +  za)  sin  x  sin  y  sin  z  + 

x  cos  x  sin  y  sin  z  +  y  sin  x  cos  y  sin  z  + 
z  sin  x  sin  y  cos  z  +  yz  sin  x  cos  y  cos  z  + 

xz  cos  x  sin  y  cos  z  +  xy  cos  x  cos  y  sin  z.  (21) 
Closed  nodal  surfaces  were  discovered  to  exist  in  the  function 

9  -  [w^5  4  CW^J]  cos  lot  (22) 

when  C  »  (18.75^*  )~X  *  .0051*0,  and  andW^^  are  given  by  Equations  (7) 
and  (21),  respectively,  with  a  «  b  -  d  ®  1  and  o  ■  B  *  6  ■  0. 

These  surfaces  are  formed  by  the  three  planes  which  contain  the  coordinate 
axes  and  by  curved  surfaces  which  meet  these  planes  at  right  angles.  The 
curved  surfaces  are  those  non-zero  points  (x^,  y^,  z^)  which  make  the 

bracket  of  Equation  (22)  equal  to  zero.  The  intersection  of  curved  surfaces 
with  the  planes  was  found  ty  evaluating  points  (x^,  y^)  for  z  =  e  with  e 


U* 


allowed  to  approach  zero  closely  enough  to  justify  the  approximations 
sin  c  ■  c  and  cos  e  *  1,  The  nodal  pattern  in  the  x-y  plane  given  by 
z  *  e  is  presented  in  Figure  5(a)*  The  pattern  which  is  Figure  5(b) 
lies  in  the  plane  z  ■  rr  and  is  repeated  for  z  ■  ±im  with  m  /  0# 


Figure  5.  Nodal  Patterns  in  z -Planes  Passed  Through  Nodal  Surfaces. 

(a)  Z  ■  c  with  e  -*  Oj  (b)  z  ■  inm  with  m  /  0.  Function 
is  +  ,oo5Uow^)  .  o. 

Since  and  are  symmetrical  to  cyclic  rotation  of  coordinates,  the 
patterns  of  Figure  5  are  equally  applicable  with  coordinates  cycled.  Based 
upon  this  property,  a  three-dimensional  model  has  been  constructed  of  these 
patterns  as  an  aid  to  visualization  of  the  volumes  enclosed  by  the  nodal 
surfaces.  Three  photographs  of  this  model  are  presented.  Figure  6(a)  rep¬ 
resents  the  anallest  enclosed  volume.  Figure  6(b)  represents  the  volume  of 
intermediate  size  enclosing  the  first.  Figure  6(c)  allows  one  to  visualize 
all  three  nodal  surfaces  and  their  enclosed  volumes.  Other  surfaces,  appar¬ 
ent  from  Figure  5*  are  not  closed  and  only  their  intersections  with  the 
principal  planes  are  shown. 

SIGNIFICANCE  OF  DISCOVERY  OF  CLOSED  VOLUMES 


In  the  theoxy  of  acoustics  the  boundary  condition  of  an  acoustic  cavity  is 
a  specified  wall  impedance,  given  by  rms  sound  pressure  in  ratio  to  rms 
particle  velocity  normal  to  the  wall.  In  symbols  this  is 

z  -  PnM-.,  (23) 

(  n)rms 

where  p  is  a  scalar  and  vR  is  one  component  of  the  vector  particle  velocity* 
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Figure  6.  Three-Dimensional  Model  Constructed  from  Patterns  of  Figure 
(a)  Innermost  Surface  Represented 
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j.aL!v 


- - 


Figure  6(b) 0  Inner  Two  Surfaces  Represented. 
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Most  frequently  the  assumed  condition  is  that  Z  is  infinite  through  the 
vanishing  of  the  normal  component  of  particle  velocity  upon  the  wall*  The 
other  two  of  the  three  mutually  orthogonal  components  of  particle  velocity 
are  not  necessarily  zero  at  the  wall  and  are  not  determined  by  that  portion 
of  wall. 


It  has  not  proved  possible  during  the  time  of  this  contract  to  discover  a 
nonseparable  velocity  potential  from  which  to  obtain  v  =  0  on  a  closed 
surface.  There  is,  however,  a  mathematical  meaning  for  the  alternate  con¬ 
dition  of  zero  impedance  upon  the  walls.  It  is  that  tne  pressure  there  is 
zero.  Tnus  we  can  say  from  Equation  (22)  that  cp  =  p,  or  that  cp  is  itself 
a  velocity  potential  from  which  we  may  derive  p  according  to 


dcp 

P  =  “  P3t  » 


(2U) 


where  p  is  tiie  mean  mass  density  of  the  medium  within  the  cavity. 

One  may  compute  frequency  of  resonance  of  the  volumes  of  Figure  6  by  re¬ 
ferring  those  volumes  to  the  cubes  which  would  enclose  them.  The  method 
is  essentially  that  described  for  membrane  patterns.  If  the  frequencies 
of  the  reference  cube  and  of  the  distorted  cube  be  u>r  and  u>c,  respectively, 

the  ratio  M  -  ‘“A0.-  is  as  follows:  for  innermost  volume,  M  *  0.983;  inner 

V  X 

two  volumes,  M  *  2.01*2;  all  three  volumes,  M  ■  3.103.  In  each  case,  the 
reference  cube  is  one  whose  side  length  is  equal  to  the  side  length  of  the 
desired  volume. 


While  enclosing  a  volume  with  walls  of  zero  impedance  is  mathematically 
justified,  it  is  nonetheless  physically  impossible.  To  produce  a  pressure- 
release  boundary  would  require  an  enclosing  wall  of  no  mass  or  stiffness, 
yet  a  wall  which  would  be  capable  of  confining  a  liquid  or  gas  without 
escape  of  matter. 

Therefore  it  is  acknowledged  that  the  significance  of  the  discovery  of 
closed  volumes  is  its  contribution  to  progress  in  understanding  the  problem 
of  applying  nonseparables  to  acoustic  cavities. 
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APPLICATION  TO  ELECTROMAGNETIC  CAVITIES 


An  electromagnetic  field  is  characterized  by  field  vectors  E  and  H.  Within 
a  lossless  cavity  these  vectors  must  individually  satisfy  the  vector  Helm¬ 
holtz  wave  equation  which  is  Equation  (1),  ana  on  the  walls  these  vectors 
must  satisfy  specif'.-  conditions.  The  conditions  consiaerea  in  this  work 
were  that  the  tangential  components  of  the  electric  field  £  and  the  normal 
component  of  the  magnetic  field  H  be  zero. 


Auxiliary  functions  known  as  the  electric  scalar  potential  and  the  magnetic 
vector  potential  are  related  to  E  and  H  everywhere  through  the  relationships 
(reference  6,  page  86) 

E  ■  -  yep  -  (25) 

H-jjVXA,  *  (26) 

and 

V  .  A  +  ue||  «  0,  (27) 


where 

tp  is  the  electric  potential 
X  is  the  magnetic  potential 

are  permittivity  and  permeability,  respectively,  within  the  cavity. 


These  potentials  satisfy  the  scalar  and  vector  Helmholtz  wave  equations, 
respectively,  wherein 


»«* U« . 


(28) 


Equation  (27)  relates  the  electric  potential  and  all  the  components  of  A. 

In  an  exercise  detailed  in  Appendix  VI,  the  electric  potential  was  equated 
to  a  nonseparable,  after  which  the  three  terms  of  the  divergence  of  A  were 
assumed  to  be  equal  to  one  another.  Components  of  A  were  then  found  by 
integration,  whereupon  field  vectors  E  and  H  follow  from  Equations  (25)  and 
(26). 


In  a  specific  case  the  electric  potential  was  chosen  to  be 


9  ■  2aVQ(x  cos  ax  sin  ay  -  y  sin  ax  cos  ay)  sin  az  eirt,  (29) 

where  aV  has  the  dimensions  of  volts  per  meter.  This  leads  to  E  and  H 
o 

components  which  are 
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where 


E„  ■  -  IjaV  cos  ax  sin  ay  sin  az  e1  . 
x  o 


E„  -  UaV  sin  ax  cos  ay  sin  az  e1  . 
y  o  3 


K  - 


u  iljk  w  <  iwt 

H  «=  -  V  sin  ax  cos  ay  cos  az  e  , 
x  3<4i  o 


i»  iiik  .  itvt 

H  ■  -  rv  V  cos  ax  sin  ay  cos  az  e  > 

y  jmjmT  oj 


i8k  iwt 

H_  -  V  cos  ax  cos  ay  sin  az  e  , 

z  3cup,  o 


•  a  to  .  a 
k  '  5c3a  • 


These  satisfy  the  reduced  Maxwell  equations  for  the  region  within  the 
cavity,  namely, 

7  .  H  -  0, 


7  •  E  -  0,  (32] 

and  they  satisfy  the  boundary  conditions  on  a  cavity  having  sides  L,  i/  and 
/  / 

L  which  obey  the  relationship 

L  t  L  *  L  ■  n  :  n  s  n  (33) 

for  n-values  that  are  any  set  of  three  positive  integers# 


Frequency  is  given  by 


,,.a  3  /IJtt  \ 

*  "  iTlV*  * 
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STUDY  OF  ELASTIC  SOLIDS 


Investigation  of  application  to  elastic  solids  has  grappled  with  two  ex¬ 
amples  of  two  boundary  conditions  at  each  edge.  In  an  effort  to  simplify 
the  elastic  problem  as  much  as  possible,  an  example  of  plane  stress  or 
plane  strain  was  first  sought. 


Consideration  was  given  to  extensional  vibration  in  a  thin  rectangular- 
plate,  a  problem  which  is  posed  and  formulated  by  Love  (reference  1,  page 
li97).  One  shear  and  one  normal  stress  are  specified  along  every  edge. 
Study  soon  showed,  however,  that  this  problem  is  mathematically  similar  to 
the  problem  of  the  right  circular  solid  cylinder  with  symmetry  about  the 
axis  (reference  h).  Nonseparable  solutions  had  already  proved  of  no  avail 
in  solving  the  cylinder  problem,  and  a  review  confirmed  the  previous  find¬ 
ing.  The  conclusion  drawn  from  consideration  of  the  extensional  vibration 
in  thin  rectangular  plates  is  that  the  inability  of  nonseparable  solutions 
to  meet  the  required  pair  of  conditions  is  other  than  a  deficiency  of  non¬ 
separable  solutions,  since  separable  solutions  do  not  meet  them  either. 


Consideration  was  later  given  to  flexural  vibration  of  a  thin  plate.  A 
rectangular  outline  was  chosen  because  all  fruitful  experience  under  the 
contract  had  been  in  rectangular  coordinates  and  because  it  had  recently 
been  fairly  well  established  that  one  cannot  have  nonseparable  solutions 
in  polar  coordinates  r  and  fl  (see  Appendix  VII).  The  flexural  wave 
equation  is  (reference  2,  page  209) 


v*n+  3P(1 zJLlta 
Eh3  at3 


0, 


(35) 


where  the  independent  variables  are  x,  y  and  t. 
Upon  assuming  that 


T) 


(36) 


one  can  write  the  reduced  flexural  wave  equation  in  factored  form  as 


where 


(Va  -  YS)(va  ♦  V)T  -  0, 


3  E 

c  -  - -r-  , 

P(1  -  *  ) 


(37) 


and 

Y  is  the  amplitude  of  displacement  at  right  angles  to  the  plate, 
h  is  the  half  thickness  of  the  plate, 

E  is  Young's  modulus, 
p  is  the  density  of  plate  material. 
v  is  Poisson's  ratio. 
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The  second  factor  in  Equation  (37)  is  the  Helmholtz  wave  equation  for  which 
separable  and  nonseparable  solutions  are  known.  The  first  factor  is  a 
Helmholtz  wave  equation  but  for  the  negative  sign.  Solutions  of  this 
equation  are  the  same  separable  and  nonseparable  solutions  that  are  already 
known,  except  that  propagation  constants  a  and  b  are  replaced  by  ia  and  ib. 

Thus  the  total  solution  of  Equation  (37)  is 


Y 


(38) 


where 


are 


with  each  a  and  b  replaced  by  ia  and  ib. 


It  is  concluded  that,  since  the  reduced  flexural  wave  equation  is  a  scalar 
equation,  functions  of  the  form 

Y  «  W<°>  +  CW^n)  +  +  DH^n\  n  =  2  (39) 

should  be  capable  of  satisfying  the  two  conditions  at  every  part  of  the 
edge,  and  closed  nodal  patterns  should  be  derivable  in  the  same  manner  as 
for  the  membrane. 
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CONCLUSIONS 


Conclusions  from  this  program  are: 

(1)  There  are  applications  for  nonseoarable  solutions  of  the  Helmholtz 
wave  equation.  They  have  been  demonstrated  for  membranes,  for  an  acoustic 
cavity  with  pressure-release  walls,  and  for  an  electromagnetic  cavity. 

They  have  been  forecast  for  a  thin  elastic  plate  in  flexural  vibration. 

(2)  Vibration  on  or  within  some  new  shapes  can  now  be  calculated  exactly 
with  functions  formed  of  the  nonseparable  solutions  added  to  separable 
solutions.  Separable  solutions  taken  alone  predict  frequency  and  dis¬ 
tribution  of  amplitude  on  squares  and  rectangles  and  in  cubes.  The  new 
shapes  are  smooth  distortions  of  the  squares  and  cubes, 

(3)  The  applicable  mathematics  is  difficult  to  handle.  The  mathematics 
proceeds  from  assumed  function  to  resulting  shape;  if  the  user  specifies 
the  shape  and  asks  the  resonant  frequency  and  amplitude  distribution,  a 
cut-and-try  effort  must  be  mounted. 

(U)  Simplifications  and  additional  applications  await  the  effort.  Insight 
gained  with  respect  to  scalar  and  vector  fields  that  use  nonseparables 
singly  and  with  separables  reveals  no  inherent  mathematical  property  which 
would  reject  this  conclusion. 
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RECOMMENDATIONS 


It  is  recommended  that  this  work  continue  with  the  objectives  of  (1)  demon¬ 
strating  more  applications  and  (2)  simplifying  the  mathematics  involved. 

Application  to  flexural  vibration  of  plates  and  to  resonant  cavities  with 
vector  conditions  on  the  walls  is  considered  attainable.  Application  to 
new  modes  in  old  shapes  and  to  prediction  of  overtones  in  new  shapes  should 
be  given  attention. 

Simplification  will  result  from  improved  organization  of  mathematical  prop¬ 
erties  which  are  recorded  in  this  report.  More  properties  await  discovery, 
among  them  being  relationships  which  express  integrals  of  no  ns  ep  arables, 
conditions  for  orthogonality,  other  effects  of  summing  variants,  and  gen¬ 
erative  operators  in  three  dimensions.  These  new  properties  will  also  con¬ 
tribute  to  simplification. 
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APPENDIX  I,  TWO-DIMENSIONAL  SOLUTIONS  IN  EXPONENTIAL  FORM 

In  two  dimensions  an  alternative  form  of  writing  the  nonseparable  solutions 
is  through  use  of  the  imaginary  i  *=  ±\£ T. 

Let 


and 


whereupon 

rewritten 


r  =  ax  t  by 

(Uo) 

s  *=  bx  +  ay, 

(Ui) 

the  generative  operator  02q  previously  written  in  a  and  b  may  be 
in  r  and  s ,  as 


°2D  c  h^a  ”  a^b  "  s37  ^2) 

The  first  five  nonseparable  solutions,  derivable  through  successive  appli¬ 
cation  of  02Dj  are 

W(0)  =  eir 
=  iseir 

*  C(is)a  -  ir]eir 

W(3)  ■■  is[(is)^  -  3ir  -  l]e^r 

W(U)  ■  C(is)4  -  6ir(is)3  +  3(ir)2  -  U(is)a  +  irle^ 

*  is[(is)4  -  lOir(is)2  +  l5(ir)2  -  10(is)a  +  l5ir  +  l]elr.  (1|3) 
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APPENDIX  II,  VARIANTS 

In  Equation  (9)  the  nonseparable  two-dimensional  solution  of  order  zero  was 
written  in  compact  form  by  introducing  phase  constants  a  and  0  to  represent 
an  arbitrary  combination  of  the  factors  in 


W  « 


(WO 


It  is  sometimes  convenient  to  write  out  explicitly  all  possible  forms  of 
the  solution,  which  in  two  dimensions  are  four  variants.  Using  numerical 
subscripts  to  distinguish  one  variant  from  another,  we  may  write 

*=  sin  ax  sin  by 
*=  cos  ax  cos  by 
*  sin  ax  cos  by 

W<°>  =  cos  ax  sin  by.  (h$) 


The  two-dimensional  generative  operator  may  then  be  applied  to  each  of  these 
with  the  result  that  four  variants  of  first  order  are  evolved.  Application 
of  O2Q  to  each  of  these  produces  four  variants  of  second  order,  and  so  on. 

Variants  of  first,  second  and  third  orders  derived  from  Equation  (U5)  are 


bx 

( -1 ) ( bx 
bx 

( -1 ) ( bx 


cos 

ax 

sin 

by  - 

ay 

sin 

ax 

cos 

by 

sin 

ax 

cos 

by  - 

ay 

cos 

ax 

sin 

hy) 

cos 

ax 

cos 

by  + 

ay 

sin 

ax 

sin 

ty 

sin 

ax 

sin 

by  + 

ay 

cos 

ax 

cos 

by) 

(-1) 


[(bx)2  +  (ay)2]  sin  ax  sin  by  +  by  sin  ax  cos 
+  ax  cos  ax  sin  by  +  2abxy  cos  ax  cos  by 

[(bx)2  +  (ay)2]  cos  ax  cos  by  -  by  cos  ax  sin 

-  ax  sin  ax  cos  by  +  2abxy  sin  ax  sin  by 

[(bx)2  +  (ay)2]  sin  ax  cos  by  -  by  sin  ax  sin 
+  ax  cos  ax  cos  by  -  2abxy  cos  ax  sin  by 

[(bx)2  +  (ay)2]  cos  ax  sin  by  +  by  cos  ax  cos 

-  ax  sin  ax  sin  by  -  2abxy  sin  ax  cos  by, 


ty 


by 


by 


by 


(U6) 


(hi) 
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*  C3(bx)2  +  (ay)2  +  l]  ay  sin  ax  cos  by  + 

2  S 

3ab(x  -  y  )  sin  ax  sin  by  + 

9  2 

3xy(a  -  b  )  cos  ax  cos  by  - 

C3 (ay)2  +  (bx)2  +  l]  bx  cos  ax  sin  by 

-  C3(ay2)  +  (bx)2  +  l]  bx  sin  ax  cos  by  + 

2  2 

3ab(x  -  y  )  cos  ax  cos  by  + 

3xy(a  -  b  )  sin  ax  sin  by  - 

[3(bx)a  +  (ay)2  +  l]  ay  cos  ax  sin  by 

■  -  C3(ay)2  +  (bx)a  +  l]  bx  cos  ax  cos  by  ♦ 

2  2 

3ab(x  -  y  )  sin  ax  cos  by  - 

3xy(a  -  b  )  cos  ax  sin  by  - 

C3(bx)2  +  (ay)2  +  l]  ay  sin  ax  sin  by 

«  [3(ay)a  +  (bx)2  +  l]  bx  sin  ax  sin  by  + 

2  2 

3ab(x  -  y  )  cos  ax  sin  by  - 
2  2 

3xy(a  -  b  )  sin  ax  cos  by  + 

C3(bx)2  +  (ay)2  +  l]  ay  cos  as  cos  by.  (U8) 

An  interesting  property  of  the  variants  is  that  the  oddness  or  evenness  of  a 
given  variant  with  respect  to  any  of  the  coordinates  is  preserved  through 

all  orders.  For  example,  the  variant  is  odd  in  x  and  even  in  y.  It 

is  found  that  '  is  likewise  odd  in  x  and  even  in  y. 

It  is  also  interesting  to  note  that  closed  nodal  lines  have  been  found  only 
with  sums  containing  like  variants,  namely,  for 

+  cwj1*  -  0,  (U9) 

and  for  one  experimental  thrust  with  second  variants.  The  last  was  the 
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function 


(a  «=  b 


•m) + 


i 

7n 


2b 


2) 


m 


(b  -  2a  -  2) 


}  -  0.  (5°) 
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APPENDIX  III,  OTHER  TWO-DIMENSIONAL  PATTERNS 

This  section  presents  a  collection  of  nodal  patterns  in  which  closed  areas 
are  not  found.  One  therefore  must  imagine  a  membrane  of  infinite  extent 
upon  which  displacement  T|  is  distributed  in  accordance  with  the  expression 

T)  *  w^e^^ »  n  *  1,  2,  3.».  (51) 

where  j  ■  1,  2,  3*  U  designates  the  variant,  and  the  patterns  are  the  curves 
upon  which  T)  ■  0  for  all  t.  Propagation  constants  in  all  cases  but  the  last 
are  a  ■  b  ■  1  so  that  for  them  frequency  is  given  by  the  frequency  equation 


0  are  presented  in  Figure  7,  in  which  variant  is  varied 


horizontally  and  order  is  varied  vertically.  The  figure  omits  the  pattern 
for  j  ■  U,  since  it  is  that  of  j  ■  3  under  a  90*  rotation  of  axes.  The  same 
scale  is  used  throughout. 


Pattern  similarities  exist  by  row  and  by  column,  particularly  in  patterns 
of  the  same  variant.  It  may  also  be  noted  that  the  straight-line  diagonals 
vanish  in  second  order,  which  may  indicate  that  all  even  orders  do  not  con¬ 
tain  them. 


Patterns  of  + 

for  n  «  1,  and  the  lower  is  for  n  ■  2.  The  distinctive  feature  of  the  pair 
is  that  it  consists  of  straight  lines  inclined  at  +U5*  to  the  horizontal 
axis  and  cut  by  one  line  at  4i5*.  The  equations  of  the  left  and  right 
patterns  are,  respectively, 

-  -  (x  +  y)  sin  (x  -  y)  (£3) 


0  are  presented  in  Figure  8.  The  upper  pair  is 


and 


w(l)  W(D 

"3  “  % 


(x  +  y)  cos  (x  -  y). 


m 


These  patterns  are  a  particularly  good  example  of  the  value  of  changing 
variables  in  the  expressions  far  nonseparables  when  a  ■  b  ■  1.  If  we  let 


P-x  +  y,  q-x-y, 

the  variants  can  be  written 
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Figure  8.  Patterns  of  -  0 

on  an  Infinite  Membrane. 
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^  (q  sin  p  -  p  sin  q) 

^  (q  sin  p  +  p  sin  q) 

^  (p  cos  q  +  q  cos  p 

wj^  ■  ^  (p  cos  q  -  q  cos  p),  *  (55) 

W^2^  •  ^  [p  sin  p  -  q  sin  q  +  p3  cos  q  -  qa  cos  p] 

W^2^  «  ^  [qa  cos  p  +  pa  cos  q  -  p  sin  p  -  q  sin  q] 

W.^2-  -  ^  Cqa  sin  p  +  pa  sin  q  +  q  cos  q  +  p  cos  p] 

^  [q*  sin  P  -  P3  sin  q  -  q  cos  q  +  p  cos  p],  (56) 

and 

-  -  \  C(q3  +  q)  sin  p  -  (p3  +  p)  sin  q  +  3pq  (cos  p  -  cos  q)] 

^  C(q3  +  q)  sin  p  +  (p3  +  p)  sin  q  +  3pq  (cos  p  +  cos  q)] 

■  -  2  [(q3  +  q)  cos  p  +  (p3  +  p)  cos  q  -  3pq  (sin  p  +  sin  q)] 

-  -  2  C(q3  +  q)  cos  P  -  (p3  ♦  P)  cos  q  -  3t  sin  p  -  sin  q)]  .  (57) 

This  change  of  variable  has  made  easier  the  task  of  calculating  patterns; 
however,  the  coordinates  of  x  and  y  have  been  used  for  all  patterns. 

Patterns  resulting  from  sums  of  first  and  zero  order  have  been  explored  to 
a  small  extent.  Closed  areas  were  found  in  two  cases  displayed  in  Figure 
2.  Figure  9  presents  five  more  patterns  of  which  the  unner  two  are  of 
the  function 

w|°)  +  CW^1)  -  0, 

with  C  ■  (5")-1  and  C  ■  (O.Wjn)”1  on  the  left  and  right  respectively.  The 
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(a)  +  (S.OOtt)’1  -o  a  -  b  -  1 

(b)  W<0)  4  (O.hliTT)’1  -  0  a  -  b  -  1 

(c)  W^0)  +  W^0)  +  -  W^)  -  0  a  -  b  -  1 


(d)  «■  (2TT)"1  -  0 


a  ■  b  •  1 


Figure  9»  Nodal  Patterns  on  Membranes  for  Various  Sums  of  First  and 
Zero  Order  Functions.  (Scales  are  the  samej  different 
portions  of  field  have  been  evaluated  in  the  five  cases.) 
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lower  three  are  of  functions  of  three  different  compositions,  one  of  which 
exhibits  closed,  egg-shaped  areas  at  some  distance  from  the  origin. 

Many  directions  for  fhture  exploration  suggest  themselves  after  study  of 
Figure  9.  One  wonders,  for  example,  if  the  combination  of  second  variants  - 
lower,  center  -  could  produce  closed  areas  with  the  appropriate  value  of  C$ 
or  again,  one  wonders  what  has  been  opened  by  the  use  of  different  propa¬ 
gation  constants  which  still  satisfy  the  frequency  equation. 
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APPENDIX  IV,  OPERATORS 


A  distinction  will  be  made  between  two  types  of  operators  encountered  in 
this  work.  One  type  is  the  generative  operator  used  to  obtain  solutions 
of  higher  order,  and  the  other  is  the  customary  partial  differential  op¬ 
erator  like  d/dx,  d/dy,  etc.  Of  the  two,  the  generative  operator  has  played 
a  greater  role  in  the  stucty-  of  nonseparable  solutions. 

Generative  Operators 

These  operators  contain  the  propagation  constants  as  variables  and,  when 
applied  to  an  nth  order  nonseparable  solution  in  rectangular  coordinates, 
yield  a  solution  of  order  n  +  1.  The  two-  and  three-dimensional  operators 
have  the  farms 

°2D  "  *3a  "  a3b 
and 

°3d-4-4-^.  <58) 


respectively, 

where 


f*=b  +  d,  g  =  a+d,  h  *=  a  -  b. 

With  the  help  of  these  operators,  a  one-to-one  correspondence  can  be  es¬ 
tablished  between  variants  of  different  orders.  If  a  variant  number  is 
assigned  to  a  zero  order  solution,  the  same  number  can  be  assigned  to  the 
highest  order  solutions  obtained  by  successive  application  of  the  operator. 
This  correspondence  is  also  preserved  when  phase  angles  or  translations  of 
coordinates  are  introduced. 

Zero,  first  and  second  order  solutions  that  are  related  through  are,  in 
compact  form,  ^ 

■  sin  u  sin  v  sin  w,  (59) 

■  fx  cos  u  sin  v  sin  w  - 

gy  sin  u  cos  v  sin  w  -  hz  sin  u  sin  v  cos  w,  (60) 

■  C(fx)a  +  (gy)a  +  (hz)3]  sin  u  sin  v  sin  w  + 

(g  +  h)  x  cos  u  sin  v  sin  v  -  2ghyz  sin  u  cos  v  cos  w  + 

(f  -  h)  y  sin  u  cos  v  sin  w  +  2fhxz  cos  u  sin  v  cos  w  + 

(f  +  g)  z  sin  u  sin  v  cos  w  2fgxy  cos  u  cos  v  sin  w,  (6l) 
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where 


u*=ax+(y,  v  «  by  +  8,  w  «=  dz  +  6. 

A  three-dimensional  nonseparable  of  first  order,  apparently  not  related  to 
those  already  described,  was  also  synthesized  from  study  of  properties  that 
a  function  must  have  to  be  a  solution  of  the  Helmholtz  wave  equation.  This 
is 

W(1)  =  x  sin  ax(b  cos  by  sin  dz  -  d  sin  by  cos  dz)  + 

y  sin  by(d  sin  ax  cos  dz  -  a  cos  ax  sin  dz)  + 

z  sin  dz(a  cos  ax  sin  by  -  b  sin  ax  cos  by).  (62) 

Mo  generative  operator  to  derive  Equation  (62)  from  Equation  (59)  with 
ry  c  B  «  6  *  0  is  yet  known. 

Differential  Operators 


It  was  found  during  the  work  that,  in  all  the  cases  tried,  the  result  of 
taking  the  partial  derivative  of  a  nonseparable  solution  with  respect  to 
x,  y*  zj  or  a  combination  of  these,  is  a  solution  of  the  wave  equation. 


Examples  are  given  below  to  illustrate  the  effects  of  differential  oper¬ 
ators  on  various  solutions.  In  two  dimensions,  the  following  relations 
hold  for  zero  order  and  first  order  solutions: 


D 

Uxwl 


DW(0) 

yi 


DxxWl0) 


-  a  w: 


3,T(0) 


DxyWl0)  "  abW2°^ 


D  ■  -  baW^^ 

yy  i  i 


j1)  .  -  aW^1)  +  bw£°) 


(1) 


-  a^1)  -  2atw(°) 


DyW^  -  bW^1^  -  aW^0) 


D  W(1) 

DyyWl 


-  b3W^^  +  2atW^0^ 


-  atWg1^  +  (b3  -  aa)W^ 

-  a3W^  -  3aabW^0)  DyyyW^1)  -  -  +  3abaW^ 

D  -  -  aabW*1^  -  (2aba  -  a3)w!^ 

xxy  1  3  '3 


(63) 
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DxyyWi1)  '  -  abi>WU1)  *  t2*2”  -  b3)wU0> 

‘  a‘wi1>  *  ‘*a3bWl°) 

D  wj1^  «  a2b2W.^  -  2ab(aa  -  b2  )VL  \ 
rxyy  1  1  v  ’  1 


D  ■  bV^  -  iiab3W^0^ 

yyyy  1  1  1 


(6U) 


For  the  examples  in  three  dimensions,  it  is  convenient  to  drop  the  variant 
notation  and  to  assign  values  to  the  phase  angles  n9  A,  and  60  The  zero- 
order  solution  from  which  solutions  of  higher  order  are  derived  is 


W 


(0) 


sin  (ax  +  a)  sin  (by  +  8)  sin  (dz  +  6). 


Unless  otherwise  indicated,  the  phase  angles  are  to  be  assumed  equal  to 
zero  in  the  following  functions : 


Dxw(0)  "  [aW'0)]  „ 

Dyw(o> .  jy°>]  „ 

D,W<°>  -  [dW<°>]  „ 


D  W^0)  -  -  aV0^ 
xx 


D  *=  -  b2w(°) 

yy 


D  =  -  dV0^ 

zz  1 


(65) 


+  (b  +  d)W^°^J  ^  “  “  a  “  2a^b  +  ^ 


v  -  "2 


DyW^ 1  ^  -  (a  +  d)W^0^]  v  DyyW^1^  B  +  2b(a  +  d)v/°) 

^  "  7 

DZW(1)  -  [dW(1^  -  (a  -  b)W(0)]  n  DzzW(1)  -  -  d2W(1)  +  2d(a  -  b)W(0\ 


(66) 


DW<2>  -  [aW^  +  2(b  +  d)^1^  -  (2a  -  b  +  d)W^0^] 

*  “  7 


D  2 ^  -  [bW^  -  2(a  +  dJW^  -  (2b  +  d  -  a)v/°)] 


8  ■ 


hO 


DzW^  -  [d1/2)  -  2(a  -  b)W^  -  (2d  +  a  4 

=  -  aV2)  -  U(b  +  d )W^ 1  ^  -  [2(b  +  d)a  - 
D  v/2^  «  -  bV2^  +  hb(a  +  d)W^^  -  [~2(a  +  d)a  - 

yy  ^ 

Dzzw(2)  -  -  dV2^  +  Ud(a  -  b)W^  -  [2(a  -  b)a  - 

Thus  in  general  the  partial  differential  operator 
operating  upon  a  nonseparable  of  order  n,  produce 
which  contains  order  n  plus  all  lesser  orders. 


TT 

7 


2a(2a  -  b  +  d)]w^0^ 


2b(2b  +  d  -  a)  tf(0) 


2d(2d  +  a  +  b)]w(0).  (6?) 

of  any  degree  will,  when 
a  sum  of  nonseparables 
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APPENDIX  V,  OTHER  THREE-DIMENSIONAL  PATTERNS 

The  function  is  derived  fran  with  two  applications  of  the  operator 

0-^  defined  in  Appendix  IV.  The  first  variant  (a  *  fl  «  fi  *  0)  with  a  =  b 

=  1  has  been  investigated  for  nodal  patterns  in  several  x-y  planes.  When 
this  is  written  in  a  manner  to  illustrate  most  clearly  the  influence  of  z 
upon  the  nodal  pattern, 

*  2  sin  z  l2(x2  +  y8  +  z  cot  z)  sin  x  sin  y  +  x  cos  x  sin  y  + 

y  sin  x  cos  y  +  i*xy  cos  x  cos  y]  *=  0.  (68) 

(2) 

From  this  it  can  be  seen  that  Wv  '  =  0  at  z  ■  0  for  all  x  and  y  and  that 
a  grid  pattern  of  sin  x  sin  y  <=  0  appears  at  z  =  tt.  Intermediate  patterns 
can  be  found  by  finding  the  (x,  y)  roots  which  make  the  square  bracket  equal 
to  zero  for  each  value  of  z  and  hence  of  z  cot  z.  The  following  table  shows 
the  substitution  to  be  made. 


VALUES  OF  Z  COT  Z  FOR  SELECTED  VALUES  OF  Z 


z 

z  cot  z 

z 

z  cot  z 

Approaching  0 

1 

±"A 

"A 

±0.8" 

-0.8tt  cot  0,2tt  ■=  -3.1*6 

±n/2 

0 

±0.85" 

-0.85"  cot  0.15"  ■  -5.21* 

±3"A 

-3"A 

±0.9" 

-0.9"  cot  0.10"  *  -8.70 

±TT 

+09 

±0.95" 

-0.95"  cot  0.05"  **  -18. 81* 

±5"A 

5"A 

±1.05" 

+1.05"  cot  0.05"  ■  +20.83 

±1.1*17" 

l 

±1.1" 

+1.1"  cot  0.10"  *  +10.61* 

±3"/2 

0 

In  approaching  the  selection  of  z  values  to  display  transition,  there  were 
few  guidelines.  As  a  consequence,  the  procedure  was  to  select  a  value, 
evaluate  a  pattern,  gauge  its  probable  interpolative  position,  and  choose 
another  value.  Therefore  the  z  values  actually  used  only  partially  cor¬ 
respond  to  those  of  the  preceding  table. 
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Nodal  patterns  which  have  been  calculated  are  presented  in  Figure  10(a)  and 

Figure  10(b).  Each  pattern  is  for  Wv  '  *=  0  in  an  x-y  plane  at  a  particular 
value  of  z  between  0  and  *,  and  the  x  and  y  values  of  each  pattern  range 
from  -3.5*  to  +3.5*.  The  x-y  plane  through  z  =  0  is  a  nodal  surface,  as 
are  the  y-z  and  z-x  planes  through  x  =  0  and  y  *=  0,  respectively.  Wavy 
nodal  surfaces  rise  wall-like  from  the  floor,  i.e.,  from  the  x-y  plane 
through  z  ■  0,  change  but  little  in  reaching  z  =  3*A>  and  then  alter 
drastically  in  finally  connecting  with  the  grid  pattern  of  z  «  tt. 

REPEATED  PATTERNS 


Consideration  of  Equation  (68)  and  the  table  of  values  of  z  cot  z  which  are 
inserted  into  Equation  (68)  shows  that  a  pattern  at  a  given  value  of  z  will 
be  repeated  for  all  other  values  of  z  which  give  the  same  value  of  z  cot  z. 
From  the  table,  the  locations  of  two  duplicate  patterns  are  shown  to  be  at 
z  approaching  0  and  at  z  *  I.I4I7*,  where  z  cot  z  ■  1,  and  at  z  «  */2  and 
z  *  3*/2,  where  z  cot  z  *=  0. 

A  general  curve  which  may  be  used  to  locate  all  identical  patterns  is  pre¬ 
sented  in  Figure  11.  This  is  a  plot  of  z  cot  z  versus  z  which  enables  one 
to  generalize  the  z-scale  through  the  use  of  the  parameter  n  «  1,  2,  3...  • 

For  example,  let  the  z-location  near  z  «  13.5  ■  U.3*  of  a  pattern  identical 

with  the  pattern  at  z  -  U.O  =  l,27n  be  desired.  In  the  organization  of 

Figure  11,  z  =  1.27*  lies  between  z  -  *  and  2tt  so  that  n  ■  2,  and  z  -  U.3n 
lies  between  z  «  U*  and  5*  so  that  n  -  $.  It  is  then  a  simple  matter  to 

read  from  Figure  11  that  the  intersection  of  z  *  U.U2*  gives  the  same  value 

of  z  cot  z  on  the  curve  of  n  «  $  that  the  intersection  of  z  ■  1.27*  gives  on 
the  curve  of  n  =  2. 


Historically,  it  was  the  pattern  of  z  «  .9919*  in  Figure  10(b)  which  led  to 
the  production  of  closed  areas  on  membranes,  for  from  it  came  recognition 
that  Equation  (68)  consisted  of  second  and  zero  order  two-dimensional  non- 
separables  (plus  a  residual  of  secondary  importance  at  lew  values  of  x  and 
y)  in  a  proportion  which  depended  upon  the  value  of  z.  Rewritten  to  em¬ 
phasize  this.  Equation  (68)  is 


z  »  const 


K  [ (w{2^  +  CvM  - 

lx  1  2D 


2  (y  sin  x  cos  y  +  x  cos  x  sin  y)J  «  0  (69) 

where  K  ■  2  sin  z  and  C  -  ^  z  cot  z  is  a  number  which  at  z  ■  .99*  has  grown 
large  enough  to  allow  the  pattern  of  to  dominate  that  of  the  sum. 
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Figure  10.  Nodal  Patterns  in  Selected 
z-Planes  for  the  Three-Dimensional. 

(2) 

Operator-Generated  Wv  7.  (Conditions 
are  a  ■  b  *  d  ■  1  and  or  ■  B  ■  6  ■  0. 
The  entire  x-y  plane  for  z  ■  0  is  a 
nodal  plane.) 


2  =  0.786377 


uu 


rjy  jH1*^  7jr'v~ 1 


Figure  11.  Generalized  Graph  of  z  cot  z  versus  z  for  Values  of  z  Between  0 
and  ±10n.  (A  horizontal  line  corresponding  to  a  given  value  of 
z  cot  z  intersects  the  curves  at  values  of  z  where  the  nodal 
patterns  are  identical.) 
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APPENDIX  VI,  VECTOR  EFFORTS 

Boundary  conditions  for  electromagnetic  cavities  are  restraints  upon  the 
electric  and  magnetic  fields  which  are  vectors.  A  common  boundary  con¬ 
dition  for  a  gas-or-liquid-fj.lled  acoustic  cavity  is  that  the  particle  ve¬ 
locity  normal  to  the  wills  be  zero,  and  particle  velocity  is  a  vector. 
Therefore  the  investigation  of  nonseparable  solutions  in  respect  to  cavi¬ 
ties  required  a  study  of  the  nonseparable s  in  relation  to  vectors. 

ACOUSTIC  CAVITIES 


This  study  began  by  treating  nonseparables  as  velocity  potentials.  Par¬ 
ticle  velocity  is  obtained  by  taking  the  gradient  of  the  potential.  For 
rigid  wall  the  requirement  is  that  a  closed  surface  exist  upon  which  the 
gradient  is  zero,  and  hence  each  component  of  the  gradient  must  simulta¬ 
neously  be  zero.  The  x-components  of  particle  velocity  for  three  cases  are 


«2D)(^3) 


o5t 


a 


[w^V,®)  “  W2D)(a/'P>]  '  2bW^)(fr/,B), 


r(D 


aw^'(or,8,6)  ..  ,  (q)  , 

vx  -  - awW(cv  ,3,6)  ♦  (b  ♦  d)W<°>(cr  ,3,6), 


(70) 

(71) 


(72) 


The  velocity  potentials  W^(ar,B)  and  W^^(fy,8,6)  are  operator-derived  from 

zero  order  in  two  and  three  dimensions,  respectively.  Application  of  the 

partial  derivative  changes  a  to  a  ■  a  -  -j.  The  propagation  constants  a, 

b,  d  are  related  through  the  frequency  equation  u>/c  ■  a  +  b  +  d,  (The 
constant  d  is  zero  in  Equation  (70).) 

The  velocity  potential  to  is  the  symmetric  function  given  in  Equation  (22). 

Nodal  patterns  of  Equations  (70)  and  (72)  were  evaluated,  the  former  under 
the  conditions  that  or  ■  3  ■  0  and  a  -  b  ■  1,  and  are  presented  as  Figures 
12(a)  and  12(b),  respectively.  In  Figure  12(a)  no  closed  areas  are  present. 
The  pattern  is  symmetrical  about  the  x  and  y  axes  but  not  about  the  ±hS* 
lines  through  the  origin.  In  Figure  12(b)  every  nodal  line  is  closed. 

This  x-y  pattern  repeats  for  z  ■  ±mrr  with  m  f  0.  At  m  -  0  the  entire  x-y 
plane  is  a  nodal  surface.  The  pattern  is  symmetrical  about  the  x  and  y 
axes  but  not  about  the  +U5*  lines  through  the  origin.  Nodal  patterns  of 
v  and  vz  at  z  -  imn  differ  from  Figure  12(b),  and  no  surface  has  been 

found  which  satisfies  the  requirement  that  the  gradient  be  zero  upon  it. 
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(a)  (b) 

Figure  12.  Nodal  Patterns  of  x-Component  of 
Velocity. 

(a)  yx  given  by  Equation  (70); 

(b)  vx  given  toy  Equation  (72)  in 
x-y  planes  given  by  z  ■  ±wtt,  m  /  0. 

ELECTROMAGNETIC  CAVITIES 


The  development  by  which  E  and  H  field  components  turned  out  to  be  separable 
solutions  when  the  electric  potential  was  nonseparable  is  given  in  general 
form. 

We  let 

„  .  W(n)ei“t,  (73) 


where  W'  is  the  operator-derived  nonseparable  of  nth  order,  and  choose  to 
examine  a  special  case,  one  in  which 

SAx  .  .  8As 

dx  dy  6z  * 
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Equation  (27)  may  then  be  written 


Therefore, 


3A  3a  SA  .  t  \  . 

*  _  y  _  Z  _  imt  w(n)  it»t 

^5y - 3z  — 3  w  e  • 


eimt  \  W<n>dx, 


(75) 

(76) 


and  similarly  for  Ay  and  Az«  (Constants  of  integration,  including  functions 

of  variables  other  than  the  variable  of  integration,  are  dropped  as  incap¬ 
able  of  satisfying  the  same  Helmholtz  wave  equation  that  w(n/  obeys,) 

Then,  by  Equation  (25),  the  electric  field  components  are 


Ex  -  -  {  tb?  ♦  t  $ w<n)dx } ei,#t 

By  -  -  {  ♦  y  $  V(n)dy  }  ei,ot 

and,  by  Equation  (26),  the  magnetic  field  components  are 


ika  (  a  e 

-  3y  J 

W(n^dz 

-  3  ' 
3z  « 

^  w(nV  ] 

feiut 

ika  {  a  C 
“lu£ J 

W^dx 

-i! 

y  W^dz  } 

loot 

e 

ika  (  &  r 
“  lou  l  7:c  J 

W^dy 

}  W^dx  } 

.  i'wt 

6  • 

(77) 


(78) 
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APPENDIX  VII,  NONSEPARABLE  SOLUTIONS  IN  POLAR  COORDINATES 

Nonseparable  solutions  were  first  found  in  circular  cylindrical  coordinates 
during  work  in  which  symmetry  about  the  axis  was  assumed.  During  this  con¬ 
tract  nonseparable  solutions  in  two  and  three  dimensions  in  rectangular  co¬ 
ordinates  have  been  explored.  Since  sane  model  configurations  would  be 
best  expressed  in  polar  coordinates,  e.g.,  thin  circular  plates,  it  was  de¬ 
cided  to  seek  nonseparable  solutions  in  r  and  0» 

The  form  of  the  separable  solution  in  r  and  0,  i.e.,  cos  (mflJJ^pr), 

suggests  that  a  nonseparable  solution  would  consist  of  a  sum  of  products  of 
functions  with  r  to  some  power  and  9  to  some  power  appearing  explicitly  in 
the  various  products. 

In  an  exploratory  effort,  the  two-dimensional  in  rectangular  coordi¬ 
nates  was  transformed  directly  into  polar  coordinates.  It  was  then  of  the 

form 

w(l)  „  tr  cos  0  cos  (ar  cos  0)  sin  (br  sin  0)  - 

ar  sin  0  sin  (ar  cos  9)  cos  (br  sin  0) 

QO 

-  br  cos  9  ^JQ(ar)  +  2  l  <**  2n(ar)  cos  2n@]  • 

n-1 

00 

L2  Z  J2n+l^  sin  ^2n  +  1^e]  " 

n«0 

00 

ar  sin  9  j^2  £  (-l)nJ2n+1(ar)  cos  (2n  ♦  1)9]  • 
n»0 

00 

[jQ(br)  +  J2n(br)  cos  2n9]*  (79) 

n-1 

This  function  was  verified  as  a  solution  of  the  wave  equation  in  polar  co¬ 
ordinates. 

However,  all  efforts  to  synthesize  a  polar  nonseparable  solution  of  two 
terms  in  analogy  to  the  first  order  in  rectangular  coordinates  have  failed. 
The  difficulty  may  stem  from  the  fact  that  the  functions  which  would  possess 
the  characteristics  suggested  in  the  second  paragraph  would  be  multivalued 
at  every  part  in  the  r,  0  plane  because  of  the  cyclical  nature  of  0. 
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